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1. Introduction 



It has long been recognized that a frequency standard 
could be based on the 282 nm transition between the 
ground 5d^°6s ^Si/2 level and the metastable 5d'^6s^ ^Dsq 
level of Hg* [I]. The lifetime of the upper level is 86(3) 
ms [2], so the ratio of the natural linewidth Ai^ to the 
transition frequency t-o is 2 X 10 '^ (Unless otherwise 
noted, all uncertainties given in this paper are standard 
uncertainties, i.e., one standard deviation estimates.) 
Doppler broadening can be avoided if the transition is 
excited with two counter-propagating photons, as origi- 
nally proposed by Bender et al. [1] and subsequently 
demonstrated by Bergquist et al. [3]. However, optical 
Stark shifts are greatly reduced if the transition is driven 
instead with a single photon by the electric-quadrupole 
interaction. In this case, Doppler broadening can be 
eliminated if the ion is confined to dimensions much 
less than the optical wavelength, as was first demon- 
strated by Bergquist et al. [4]. 



Recently, the (F = 0, Mf = 0) to {F = 2, Mf = 0) hy- 
perfine component of the "'Hg* 5d^°6s ^Sic to 
5d'^6s^ ^D5/2 single-photon transition has been observed 
with a linewidth of only 6.7 Hz by Rafac et al. [5]. A 
laser servo-locked to this transition is an extremely sta- 
ble and reproducible frequency reference. New develop- 
ments in optical frequency metrology [6, 7] may soon 
make this system practical as an atomic frequency stan- 
dard or clock. 

While the (F = 0, Mf = 0) to (F =2,Mf = 0) hyper- 
fine component has no linear Zeeman shift, it does have 
a quadratic Zeeman shift that must be accounted for. In 
addition, there is a second-order Stark shift and a shift 
due to the interaction between the electric-field gradient 
and the atomic electric-quadrupole moment. None of 
these shifts has yet been measured accurately, so it is 
useful to have calculated values, even if they are not very 
precise. Also, it is useful to know the functional form of 
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the perturbation, even if the magnitude is uncertain. For 
example, the quadrupole shift can be eliminated by aver- 
aging the transition frequency over three mutually or- 
thogonal magnetic-field orientations, independent of the 
orientation of the electric-field gradient. 



2. Methods and Notation 

The quadratic Zeeman shift can be calculated if the 
hyperfine constants and electronic and nuclear g -factors 
are known. Similarly, the quadratic Stark effect can be 
calculated from a knowledge of the electric-dipole oscil- 
lator strengths. The quadrupole shift depends on the 
atomic wavefunctions. Some of these parameters have 
been measured, such as the hyperfine constants and 
some of the oscillator strengths. There are also pub- 
lished calculations for some of the oscillator strengths. 

Here, we estimate, by the use of the Cowan atomic- 
structure codes, values for parameters for which there 
are neither measured values nor published calculations. 
The Cowan codes are based on the Hartree-Fock ap- 
proximation with some relativistic corrections [8]. The 
odd-parity configurations included in the calculation 
were 5d'°np (« = 6,7,8,9), 5rf"'5/, 5d%s6p, 5d%slp, 
5d^6s5f, and 5d^6s^6p. The even-parity configurations 
were 5d'"ns (« = 6,7,8,9,10), 5d'°nd (n = 6,7,8,9), 
5d'^6s^, 5d''6s7s, 5d''6s6d, and 5d''6p^. Recently, San- 
sonetti and Reader have made new measurements of the 
spectrum of Hg* and classified many new lines [9]. 
They also carried out a least-squares adjustment of the 
energy parameters that enter the Cowan-code calcula- 
tions in order to match the observed energy levels. We 
use these adjusted parameters in our Cowan-code calcu- 
lations. 

As one test of this method of calculation, we esti- 
mated the weakly allowed 10.7 /j,m 5J"'6/? ^Pic to 
5d''6s^ ^D3/2 electric-dipole decay rate. This decay is al- 
lowed only because of configuration mixing, since it 
requires two electrons to change orbitals. The calcula- 
tion shows the decay to be due mostly to mixing between 
the 5d"'6/7 and 5d'^6s6p configurations. The calculated 
rate is 111 s '; the measured rate is 52(16) s ' [2]. 
Another test is the electric-quadrupole decay rate of the 
5d'^6s^ ^D5/2 level to the ground level. The calculated rate 
is 12.6 s ', and the measured rate is 11.6(0.4) s '. Sim- 
ilar calculations have been carried out by Wilson [10]. 

Let Ho be the atomic Hamiltonian, exclusive of the 
hyperfine and external field effects, which are treated as 
perturbations. For convenience, we denote the eigen- 
states of Ha corresponding to the electronic levels 
5d^°6s ^Si/2 and 5d''6s^ ^'Dm having 7^ eigenvalue Mj by 
IS 1/2 Mj) and ID 5/2 Mj), respectively. 



The corresponding eigenvalues of Hq are denoted 
W(S, 1/2) and W(D, 5/2). An arbitrary eigenstate of Ho 
with eigenvalue W{y, J) and electronic angular mo- 
mentum 7 is denoted \y J Mj). Since ''"Hg* has in addi- 
tion a nuclear angular momentum /, where / = 1/2, the 
complete state designation is \y JFMp), where F is the 
total angular momentum, and Mf is the eigenvalue of F^ . 



3. Quadratic Zeeman Shift 

In order to calculate the energy shifts due to the 
hyperfine interaction and to an external magnetic field 
B = Bz, we define effective Hamiltonian operators 
H's and Ho that operate within the subspaces of hyper- 
fine sublevels associated with the electronic levels 
5(i"'6.s ^Si/2 and 5d''6s^^T)si2, respectively: 

H's = hAsI-J + g.,{S)fiB JB + g',piBl-B, (1) 

Hh = hAuI-J + g.,(D)/^B JB + g',fjUBl-B , (2) 

where As and Ad are the dipole hyperfine constants, 
gj{S) and gj(D) are the electronic g-factors, gl is the 
nuclear g -factor, h is the Planck constant, and jjlb is the 
Bohr magneton. All of the parameters entering Hs and 
Hi, are known from experiments, although a more accu- 
rate measurement of gj(D) would be useful. The 
ground-state hyperfine constant As has been measured 
in a "''Hg* microwave frequency standard to be 40 
507.347 996 841 59 (43) MHz [11]. The excited-state 
hyperfine constant Ad has been measured recently by an 
extension to the work described in Ref. [5], in which the 
difference in the frequencies of the IS 1/2 0) to ID 5/2 
2 0) and the IS 1/2 0) to ID 5/2 3 0) transition frequen- 
cies was determined to be 3Ad = 2 958.57(12) MHz 
[12], in good agreement with an earlier, less precise 
measurement by Fabry-Perot spectroscopy [13]. The 
ground-state electronic g -factor g.;(S) was measured in 
'""^Hg* by rf-optical double resonance to be 2.003 174 
5(74) [14]. The excited- state electronic g -factor gj(D) 
was measured in '""^Hg* by conventional grating spec- 
troscopy of the 398 nm 5d'%p ^P^n to 5d%s^ ^Dsn line 
to be 1.198 0(7) [15]. The difference in gj(S) or gj(D) 
between '''^Hg* and '''Hg* is estimated to be much less 
than the experimental uncertainties. The nuclear g -fac- 
tor gl is -5.422 967(9) X 10"" [16]. The measurement 
was made with neutral ground-state ''"Hg atoms, so the 
diamagnetic shielding factor will be slightly different 
from that in the ion. However, this is effect is negligible, 
since the magnitude of g! is so small compared to gjiS) 
or gAD). 
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The determination of gj(D) could be improved by 
measuring the optical-frequency difference between 
two components of the 282 nm line and the frequency of 
a ground-state microwave transition at the same mag- 
netic field. Since the uncertainty in the quadratic Zee- 
man shift is due mainly to the uncertainty in gj(D), it is 
useful to see how accurately it can be estimated theoret- 
ically. The Lande g -factor for a ^D5/2 state, including the 
correction for the anomalous magnetic moment of the 
electron, is 1.200 464. The Cowan-code calculation 
shows that the configuration mixing does not change 
this value by more than about 10 *■, i.e., 1 in the last 
place. There are several relativistic and diamagnetic cor- 
rections that modify g./(D), one of which, called the 
Breit-Margenau correction by Abragam and Van Vleck 
[17], is proportional to the electron mean kinetic energy. 
The other corrections are more difficult to calculate. 
The Cowan-code result for the mean kinetic energy of 
an electron in the 5d orbital of the 5d'^6s^ configuration 
is r= 19.32 hcR^, where Ro= is the Rydberg constant. 
Using this value, we obtain a theoretical value of gj{D), 
including the Breit-Margenau correction, of 1.199 85, 
which disagrees with the the experimental value by 
1.85 X 10 ^ which is 2.6 times the estimated experi- 
mental uncertainty of Ref. [15]. If we calculate gj(D) 
for neutral gold, which is isoelectronic to Hg*, by the 
same method, we obtain a value which differs from the 
accurately measured experimental one [18] by 
(7 ± 2) X 10"^. Thus, the error in the calculated value 
for gj(D) of '"'Hg* might be less than 1 X 10 ', but it is 
impossible to be certain of this, since there are uncalcu- 
lated terms. Measurements of the '""Hg* optical clock 
frequency at different values of the magnetic field 
should result in a better experimental value for gj(D) in 
the near future. 

For low magnetic fields (B less than 1 mT), it is 
sufficient to calculate the energy levels to second order 
in 5 . To this order in B , the energies of the hyperf ine- 
Zeeman sublevels for the ground electronic level are 



W(S, 1/2,0, 0,B) = W(S, 111) 



[g/(S) - g',fti.lB' 



3hAs 



4Ms 



W(S, 1/2, 1, 0,B) = WiS, 1/2) H- 



(3) 



Ms 



[gj(S) - glffJjJB' 
4ftAs 



(4) 



W(S, 1/2, 1, ±l,B) = W(S, 1/2) H- 



[gj(S) + gnf^BB 



For the 5d''6s^ ^Ds/t level we have 



Ms 



W(D, 5/2, 2,0, B) = W(D, 5/2) - 

[gj{D) - glftiJB' 
12/iAd 

WiB, 5/2, 2, ±1, S) = W(D, 5/2) 



ThAo 



IhAo 



[7g,(D) - g;]MBfi 2[g,(D) - giff^JB' 



27/zAd 



W(D, 5/2, 2, ±2, S) = W(D, 5/2) 



7/jAd 



[7g,(D) - g;]MBfi 5[g.,(D) - giff^jB' 



3 



108Md 



W(D, 5/2, 3, 0, B) = W{D, 5/2) + 



[g,(D) - glffJjJB' 



5Md 



12Md 



W(D, 5/2, 3, ±1, B) = W(D, 5/2) + 



5hAo 



27Md 



W(D, 5/2, 3, ±2, B) = W(D, 5/2) + 



5hAo 



W(D, 5/2, 3, ±3, B) = W(D, 5/2) + 



[5gy(D)-Hg;]/XBB 

2 



5hAo 



(5) 



(6) 



(7) 



, (8) 



(9) 



[5g,(D) + g;]/^BB , 2[g,(D) - g'lffJilB' ^jQ^ 



[5g,(D) + g;]/^BB 5[g,(D) - giflJiiB' 

3 108/2Ad ' ^ ' 



(12) 



Here, W{y, J, F, Mf, B) denotes the energy of the state 
lyJFMf), including the effects of the hyperfine interac- 
tion and the magnetic field. 
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At a value of B of 0. 1 mT, the quadratic shift of the 
IS 1/2 0> to ID 5/2 2 0) transition (optical clock transi- 
tion) is —189.25(28) Hz, where the uncertainty stems 
mainly from the uncertainty in the experimental value 
of gj(D)- In practice, the error may be less than this if 
the magnetic field is determined from the Zeeman split- 
tings within the ID 5/2 F Mp) sublevels. The reason is 
that an error in gj{D) leads to an error in the value of B 
inferred from the Zeeman splittings, which partly com- 
pensates for the gj(D) error. If instead we use the calcu- 
lated value of gj{D), the quadratic shift for B = 0.1 mT 
is —189.98 Hz, where the uncertainty is difficult to 
estimate. 



4. Quadratic Stark Shift 

The theory of the quadratic Stark shift in free atoms 
has been described in detail by Angel and Sandars [19]. 
The Stark Hamiltonian is 

HE=-iJiE, (13) 

where /it is the electric-dipole moment operator. 



/t= -e 



Sn, 



(14) 



ftscalaXy, /) ■ 






l(y7ll^i-">lly'7')P 



3{2J+\)^,W{y\J')-W{y,J) 



(16) 



,r(y, /) = Siret 



107(27- 1) 



,3(27h-3)(7h-1)(27h-1). 



y j_A\ 1 2] l(y7ll/x"'lly7')P 

The summations are over all levels other than Iy7). 
Equations (16) and (17) can be rewritten in terms of the 
oscillator strengths /y .;,,,;■ : 



ascalar(y, 7) = 2j 



fyJ.' 



y'r 



,r(y, 7) = 



;.[W(y',7')- W(y,7)]- 



A-ne^e^h^ [307(27 - 1)(27h- 1) 



,(18) 



(27-h3)(7-h 1) 



x2(-iH',!?]fw.^ 



y.r 



yJ.y'J' 



J 7 J'hw{y',J') - W(y,J)f 



(19) 



where nic is the electron mass. The tensor polarizability 
is zero for levels with 7 < 1, such as the Hg* 5d"^6s ^Siq 
level. 

For an atom with nonzero nuclear spin /, the 
quadratic Stark shift of the state lyJFMf) is 



and E is the applied external electric field. In Eq. (14), 
r, is the position operator of the ith electron, measured 
relative to the nucleus, and the summation is over all 
electrons. 

First consider an atom with zero nuclear spin, such as 
i98]_[g+ 'Yo second order in the electric field, the Stark 
shifts of the set of sublevels ly7Mj) depend on two 
parameters, ascaiaXy, 7) and atcnsor(y, 7), called the 
scalar and tensor polarizabilities. In principle, when 
both magnetic and electric fields are present but are not 
parallel, the energy levels are obtained by simulta- 
neously diagonalizing the hyperfine, Zeeman, and Stark 
Hamiltonians. In practice, the Zeeman shifts are nor- 
mally much larger than the Stark shifts, so that He does 
not affect the diagonalization. In that case, the energy 
shift of the state \y JMj) due to He is 



^WE{y, 7, Me, E)=-\ a,ea,ar(y, 7, F ) E^ 



, , ^J3M^-F(F+1)] 



,r(y, 7, F) ■ 



F(2F - 1) 



(3E^ - E^). (20) 



We make the approximation that hyperfine interaction 
does not modify the electronic part of the atomic wave- 
functions (the 77 -coupling approximation of Angel and 
Sandars [19]). This approximation is adequate for the 
present purpose, which is to evaluate the Stark shift of 
the '""Hg* optical clock transition. Obtaining the differ- 
ential Stark shift between the hyperfine levels of the 
ground state, which is significant for the '^'Hg* mi- 
crowave frequency standard [11], requires going to a 
higher order of perturbation theory [20]. In the /7-cou- 
pling approximation [19], 



^WE{y, 7, M,, E)=-\ a,calar(y, J)E^ 



"tcnsoXy, 7) 



[3M| -7(7-H 1)] 
7(27- 1) 



(iEj^ - E^). (15) 



Treating He by second-order perturbation theory leads 
to the following expressions for the polarizabilities [19]: 



ascalar(y, 7, F) = ttscalaXy, 7), 

a,c„»r(y,7,F) = (-l)'*'*^ 

F(2F - \){2F -H 1)(27 -H 3)(27 -H 1)(7 -H 1)" 
(2F-h3)(F-h 1)7(27- 1) 



(21) 



\F J I 

\ J T7 •") ("tensor 



{J, J). 



(22) 
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Equations (18) and (19) were used to evaluate the 
polarizabilities for the Hg* 5d^°6s ^Sm and 5d^6s^ ^D^n 
levels. For the calculation of ascaiar(S, 1/2), the oscillator 
strengths for all electric-dipole transitions connecting 
the 5d^%s configuration to the 5d^°np (n = 6,7,8) and 
5d^6s6p configurations were included. These were 
taken from the theoretical work of Brage et al. [21]. The 
final result is a,caiar(S, l/2)/(4T7eo) = 2.41 X 10"'* cm', 
which compares very well with the value of 
2.22 X 10"^'' cm"* obtained by Henderson et al. from a 
combination of experimental and calculated oscillator 
strengths [22]. For the calculations of ascaiar(D, 5/2) and 
atcnsor(D, 5/2), the oscillator strengths for electric-dipole 
transitions to the 5d'°np {n = 6,7,8), 5d'°5f, and 
5d'^6s6p configurations were taken from Brage et al. 
[21]. The oscillator strengths for electric-dipole transi- 
tions to the 5d^6s7p and 5d'^6s^6p configurations were 
taken from the Cowan-code calculations. The results 
were ascaiar(D, 5/2)/(4TTeo) = 3.77 X 10"'" cm' and 
a,e„sor(D, 5/2)/(4TTeo) = -0.263 X lO"''* cml Evaluating 
Eq. (22) for F = 2 and F = 3 in the 5d%s^ 'D5/2 level, we 
obtain a,e„sor(D, 5/2, 2) = |ate„sor(D, 5/2) and a,e„sor(D, 
5/2, 3) = a,e„sor(D, 5/2). 

The tensor polarizability is much smaller than the 
scalar polarizabilities and in any case does not con- 
tribute if the external electric field is isotropic, as is the 
case for the blackbody radiation field. The net shift of 
the optical clock transition due to the scalar polarizabil- 
ities is ^[ascaiar(S, 1/2) — a,caiar(D, 5/2)]£''. lu frcquency 
units, the shift is —1.14 X 10"' E^ Hz, where E is ex- 
pressed in V/cm. The error in the coefficient is difficult 
to estimate, particularly since it is a difference of two 
quantities of about the same size. However, the total 
shifts are small for typical experimental conditions. If 
the electric field is time-dependent, as for the blackbody 
field, the mean-squared value (£') is taken. At a temper- 
ature of 300 K, the shift of the optical clock transition 
due to the blackbody electric field is —0.079 Hz. The 
mean-squared blackbody field is proportional to the 
fourth power of the temperature. For a single, laser- 
cooled ion in a Paul trap, the mean-squared trapping 
electric fields can be made small enough that the Stark 
shifts are not likely to be observable [23]. 



5. Electric Quadrupole Shift 



The interaction of the atomic quadrupole moment 
with external electric-field gradients, for example those 
generated by the electrodes of an ion trap, is analogous 
to the interaction of a nuclear quadrupole moment with 
the electric field gradients due to the atomic electrons. 
Hence, we can adapt the treatment used for the electric- 
quadrupole hyperfine interaction of an atom [26]. The 
Hamiltonian describing the interaction of external elec- 
tric-field gradients with the atomic quadrupole moment 



Hq = V£«'-0<'' = 2(-l)'V£f 6>i'^, 



(23) 



,=2 



where V/i'" is a tensor describing the gradients of the 
external electric field at the position of the atom, and 
0*'* is the electric-quadrupole operator for the atom. 

Following Ref. [26], we define the components of 
V£"' as 



2 dz 



VE9} = 



Ve dE^ 

6 dz 



V6 



d^E, 



VE?l=-^d._E._ 



(24) 
(25) 
(26) 



where E± = E^ ± \Ey and 5± = |^ ± ij;. 

The operator components 0,f' are defined in terms of 
the electronic coordinate operators as 



&^'^=-lJ.Ozf-rf), 



m = -e^l'-j,zMj±m\ 



m=-e\'^J.(x,±iy,f 



(27) 
(28) 
(29) 



where the sums are taken over all the electrons. The 
quadrupole moment @iy, J) of an atomic level ly /) is 
defined by the diagonal matrix element in the state with 
maximum Mj : 



The atomic quadrupole moment is due to a departure 
of the electronic charge distribution of an atom from 
spherical symmetry. Atomic quadrupole moments were 
first measured by the shift in energy levels due to an 
applied electric-field gradient in atomic-beam reso- 
nance experiments [24, 25]. 



@{y,J) = {yJJ\@lP\yJJ) 



(30) 



This is the definition used by Angel et al. [24]. 

In order to simplify the form of VE*'*, we make a 
principal-axis transformation as in Ref. [27]. That is, we 
express the electric potential in the neighborhood of the 
atom as 
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^{x', y\ z') = A [(x'^ + y'^ - 2z") + e{x'^ - y'^)]. (31) 

The principal-axis (primed) frame (x', y', z') is the one 
in which (p has the simple form of Eq. (31), while the 
laboratory (unprimed) frame {x,y,z) is the one in 
which the magnetic field is oriented along the z axis. 

The tensor components of V/i*^' in the principal-axis 
frame are obtained by taking derivatives of <P(x' , y' , z'): 



V£^^''=-2A, 



VMY = 0, 



(32) 
(33) 



VSlf = A^eA. 



(34) 



In the principal-axis frame, Hq has the simple form 



2A@P' + xl^ eA[ 0?'' + ©if ). (35) 



As long as the energy shifts due to Hq are small 
relative to the Zeeman shifts, which is the usual case in 
practice, Hq can be treated as a perturbation. In that 
case, it is necessary only to evaluate the matrix elements 
of Hq that are diagonal in the basis of states lyJFMp), 
where F is the total atomic angular momentum, includ- 
ing nuclear spin / , and Mp is the eigenvalue of F-^ with 
respect to the laboratory (not principal-axis) frame. Let 
a> denote the set of Euler angles {a, 13, y} that takes the 
principal-axis frame to the laboratory frame. To be ex- 
plicit, starting from the principal-axis frame, we rotate 
the coordinate system about the z axis by a , then about 
the new y axis by fi , and then about the new z axis by 
y so that the rotated coordinate system coincides with 
the laboratory coordinate system. We can set 7=0, 
since the final rotation about the laboratory z axis, 
which is parallel to B , has no effect. The states I y JFm ) ' 
defined in the principal-axis frame and the states 
\y JFjx) defined in the laboratory frame are related by 



\y JFmY = ^ D'^,{oy)\y JFiJi), 



(36) 



where D^f„'(w) is a rotation matrix element defined in 
the passive representation [28, 29]. The inverse relation 
is 



\yJFtJi) = 2 D';^'(M)\yJFmy. 



(37) 



In order to evaluate the diagonal matrix elements of 
Hq in the laboratory frame, it is necessary to evaluate 
matrix elements of the operators #^^' , defined in the 



principal-axis frame. These matrix elements are of the 
form 

(yJFfjLie^'^lyJFfi) 
= ^D<^Aio)D'^'iM) '{yJFm'\&f\yJFm)\ (38) 



:(y/FII(9<2>ll77F)2(-l)''""' 



XI '^ , ^ ^)d<^'.(w)£><T(w), 
-m q m/ 



■.(-lY-t^-iiyJFm^^MyJF) 



(39) 



F 2 F 

-m' q m 



x2( _' , " ' )D<^J(w)Z)i''^-„(w), 



(40) 



{-\Y-^--1{yJF\\e^'■^\yJF) 2 {2K+\) 

Kmm'nn' 

-m q m/\fji — jjl n/\m —m n j 

(41) 

= {-\i'->'-\y]F\\@'''^\yJF)[^l 2 ^^^p).^^)^ 

(42) 

where Eq. (39) follows from the Wigner-Eckart theo- 
rem, and Eqs. (40), (41), and (42) follow from Eqs. 
(4.2.7), (4.3.2), and (3.7.8) of Ref. [28], respectively. 
The required rotation matrix elements are, from Eq. 
(4.1.25) of Ref. [28] (with correction of a typographical 
error). 



D\^\m) = i(3 cos^ /3 - 1), 



(43) 



Di5?2(w) = A /I sin^ /3(cos 2a + i sin 2a). (44) 



The 3-j symbol in Eq. (42) is 



F 2 F 
-jj. iJ, 



(-1) 



F-IJ.. 



2[3fji'-F(F+l)] 



[{2F + 3)(2F + 2){2F + 1)2F(2F - 1)]"' ' 

(45) 
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The diagonal matrix elements of Hq in the laboratory 
frame are 

{y JFMfIHqIj JFMf) 

_ -2[3M| - F(F+ l)]A(yyFII6)"'lly/F) 
~ l(2F + 3)(2F + 2)(2F + l)2F(2F- 1)]'" 

X [(3 cos^ )8 - 1) - e sin^ ;S(cos^ a - sin^ a)]. (46) 

It is simple to show, by directly integrating the angu- 
lar factor in square brackets in Eq. (46), that the average 
value of the diagonal matrix elements of Hq, taken over 
all possible orientations of the laboratory frame with 
respect to the principal-axis frame, is zero. This also 
follows directly from the fact that the quantity in square 
brackets is a linear combination of spherical harmonics. 
It is less obvious that the average, taken over any three 
mutually perpendicular orientations of the laboratory z 
quantization axis, is also zero. This result is proven in 
Appendix A. This provides a method for eliminating the 
quadrupole shift from the observed transition frequency. 
The magnetic field must be oriented in three mutually 
perpendicular directions with respect to the trap elec- 
trodes, which are the source of the external quadrupole 
field, but with the same magnitude of the magnetic field. 
The average of the transition frequencies taken under 
these three conditions does not contain the quadrupole 
shift. 

The reduced matrix element in Eq. (46) is, in the 
//-coupling approximation, 

(y(/7)FII6)<"lly(/7)F) 



= {-iy*'*''(2F+ 1) 



J 2 J 
F I F 



J 2 J 
-J J. 



@(y, J), (47) 



where / is included in the state notation in order to 
specify the order of coupling of / and J. For the partic- 
ular case of the '''Hg* 5d'^6s^ ^Gsn level, the reduced 
matrix elements are 

(D 5/2 2II0<^>IID 5/2 2) = 2 J^ 0(D, 5/2), (48) 



(D 5/2 3II0<">IID 5/2 3) = 2yy 0(D, 5/2), (49) 

Since the Cowan-code calculation shows that there is 
very little configuration mixing in the ''''Hg* 5d'^6s^ ^Dsq 
level, @{D, 5/2) can be reduced to a matrix element 
involving only the 5d orbital: 



@(D,5/2) = ^{5d%,2,mj 



5/2\3z^ - r^\5d^ds,2, mj = 5/2), 



= I {5d, m, = 2l3z' - r^\5d, m, = 2), 



(50) 
(51) 

e \j^ {5d, m, = 2|}'2.o(6', (^)l5fif , m, = 2), (52) 



■■e-J^{5d\r^\5d) 



X I I ¥2.2(0, </))F2.o(0, (t>)Y2,2(9,(l})sm 0d0d(t>, (53) 
^0 •'0 



■■5e{5d\r^\5d) 



2e 
1 



2 2 2\/2 2 2 



2 2/VO 



{5d\A5d). 



(54) 



(55) 



The apparent sign reversal in Eq. (50) relative to Eqs. 
(27) and (30) is due to the fact that the quadrupole 
moment is due to a single hole in the otherwise filled 5d 
shell rather than to a single electron . According to the 
Cowan-code calculation, 

{5d\r'-\5d) = 2.324 al = 6.509 X IQ-'^ cm^ (56) 

where do is the Bohr radius. 

Since the quadrupole shifts are zero in the 5d'^°6s ^Sm 
level, the quadrupole shift of the '""Hg* optical clock 
transition is due entirely to the shift of the ID 5/2 2 0) 
state, and is given by 

(D 5/2 2 OIHqID 5/2 2 0> 

= fA6>(D, 5/2)[(3cos'/3- 1) 

— e sin^ j8(cos^ a — sin^ a)], (57) 
= - i-,Ae{5d\r^\5d)l(3 cos' (3-1) 

— € sin' ;S(cos' a — sin' a)], (58) 
= - 3.6 X lO^'M [(3 cos' [3- I) 

- € sin' j8(cos' a - sin' a)]Hz, (59) 

where A is expressed in units of V/cm'. Thus, for typical 
values A ~ 10^ V/cm' and \e\^ 1, the quadrupole shift is 
on the order of 1 Hz. 
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6. Appendix A. Angular Averaging of the 
Quadrupole Shift 

For the purpose of describing the quadrupole shift, 
the orientation of the laboratory (quantization) axis with 
respect to the principal-axis frame is defined by the 
angles /3 and a . In the principal-axis coordinate system, 
a unit vector along the laboratory z axis is defined in 
terms of /3 and a by 



z = (sin /3 cos a, sin ;S sin a, cos ;S). 



(60) 



We wish to show that the angular dependence of the 
quadrupole shift is such that the diagonal matrix ele- 
ments given by Eq. (46) average to zero, for z along any 
three mutually perpendicular directions. 

An arbitrary set of three mutually perpendicular unit 
vectors ei, ej, and e^, can be parameterized by the set of 
angles 9,4'^ and ip in the following way: 



Ci = (sin cos (f), sin sin 4>, cos B), 



(61) 



62 = (cos (f> COS COS il/— sin (f> sin t//, sin cos 6 cos ij/ 



+ cos cf) sin ip, —sin cos i//). 



(62) 



e3 = (— cos<^cos0sint//— sin <^ cos i//,— sin (^ cos 0smip 

+ cos (f> COS ij/, sin sin t//). (63) 

It can be verified by direct computation that c.-e^ = Sy. 
The quadrupole shift can be evaluated for each of 
these three unit vectors substituted for z [Eq. (60)] and 
the average taken. First consider the average of the quan- 
tity (3 cos^ 13—1) that appears in Eq. (46): We use the 
fact that cos /3 is the third component of f , so the 
average is: 

(3 cos^ j3 — I) = cos^(?H- sin^0 cos^(/f H- sin^0 sin^)// — 1, 

(64) 

(65) 

(66) 



= cos^0H- sin^e - 1, 



^0, 



for arbitrary 0, 4>, and (//. Similarly, the average of 
the other angle-dependent term in Eq. (46), 
sin^ /3(cos^ a — sin^ a), is calculated by making use of 
the fact that sin /3 cos a is the first component of f , and 
sin ;S sin a is the second: 

(sin^ j8(cos^ a — sin^ a)) 

= |[sin^ cos^ 4> ~ sin^ sin^ </> 



-I- (cos 4) COS COS ip — sin 4> sin ifsf 
— (sin (j) cos cos ij/ + cos (f) sin ij/f 
+ (cos <(> cos sin il/+ sin <^ cos if/y 
(sin (^ cos sin (// — cos <p cos i/;)^], 
= 0, 



(67) 
(68) 



for arbitrary 0, 4>, and i//. Hence, the matrix elements of 
Hq given by Eq. (46) average to zero for any three 
mutually perpendicular orientations of the laboratory 
quantization axis. 
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